QUANTIZED GEODESIC FLOW AND SUPERCONDUCTIVITY OF PLASMA IN 

FIREBALLS 



One looks for superconductivity at room temperature. But the superconductivity is 
a purely quantum phenomenon and it is worthwhile to search it in situations in which 
the quantum nature of substances is most clearly defined. These are, naturally, the low 
temperature, and besides the state of plasma for comparatively high temperature. I 
belive that plasma can be a superconductor and the verification of this suggestion is the 

£NJ ' phenomenon of fireball. 

In the work "On the nature of fireball" [13] P.L.Kapitza calculated that the whole 
. energy which can be contained in a ball with the size of a standard fireball (even for the 
case of full ionization) would irradiated during the time of order 0.01 sec. But the life 
| time of fireballs is several minutes, that has increased the theoretical value ten thousend 
times. Because of this P.L.Kapitza was forced to propose the hypothesis of constant 
■ feeding fireballs by the energy of external magnetic field. But the succeeding experiments 

Q ■ of generation such fields by Super- High- Frequency vibrations (SHF) to derive fireballs 
have been unsuccessful. There were tens of efforts to find a physical mechanism which 
could explain such a long life time of fireballs. It is generally agreed that all of them 
collapsed. It seems to me that there is the only way out of the difficulty. To accept that 
the groundwork of the phenomenon is the superconductivity of plasma. 

,— i ■ B.B.Kadomzev in the book "Collective phenomena in plasma" [12] considers oscil- 

lations in plasma with negative energy which is indirect evidence for a gap in energy 
spectrum and hence for the possibility of superconductivity. In the paper of collaborators 
of Lawrens Livermore National Laboratory "Temperature Measurements of Shock Com- 
preesed Liquid Deuterium up to 230 GPa" [17] it was described that at high pressures 
and temperatures deuterium behaves like a degenerate Fermi-liquid metal: its compress- 
ibility and conductivity abruptly increased, it appears a highly reflective state which is 
caracteristic of a liquid metal. Such behaviour also points out to a possibility of supercon- 
ductivity. An added reason for superconductivity is the Meisner-Ochsenfeld effect. The 
fact that magnetic fields penetrate into electronic plasma only to the depth of A = c/uj 
shows strong diamagnetic properties of plasma. Here c is the light velocity; u is the 
Langmure frequency (the frequency of "electronic sound"); u = ^J4:7ie 2 n/m , where n is 
the electronic density, m is the mass of electron, e is its charge [12]. The value A = c/uo 
plays the part of Debye wavelength. In the paper of Arsenjev [1] there is another indirect 
evidence in favour of superconductivity. He indicates that collisions in plasma leads to 
very small changes of momentum. 

Let us try to describe a qualitative picture of the phenomenon under the conjecture of 
superconductivity of plasma. Let us imagine ionized plasma that generates a positively 
charged ball and a spherical cover around this ball generated by negatively charged elec- 
tronic plasma. Suppose at last that in this cover it is realized superconductivity. Electrons 
with momentum which tangent to the sphere are deflected by the positive charge of the 
central part of the ball remaining in the spherical cover. The superconductivity means 
that the energy of electrons do not looses by collisions. Any continuous vector field on a 
sphere has singularities. To avoid it, suppose that in our case something like a geodesic 
flow on the sphere takes place. Preservation of the ionization state means that one deals 
with the cold plasma. This plasma does not failed due to superconductivity. It explaines 
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such a long life time of fireballs. A ball of plasma is hanged in the air and in this ball 
takes place the current. It is unusual current since it does not have a fixed direction at 
any point and at any moment. It is spreaded in all possible tangent directions. Inside the 
ball there is positive charge which directs and maintains this current. Such phenomenon 
is possible at the state of superconductivity only — a current that does not have a direc- 
tion. Electronic liquid flows indeed, but it flows simultaneouly in all directions. Hence 
the fireball is a physical realization of the mental idea of Vlasov (which was used by him 
as the basis for his equation) [8], the idea of probability distribution of electrons relative 
to its positions and velocities. 

How to organize such current at the low temperature under the usual superconductiv- 
ity? It would be interesting to take a negatively charged supercoductor of the spherical 
form and to put inside it a sufficiently large positive charge. 

Thus the fireball, having the same density as the air and being of no different from it 
except for its state, can freely float in the air which is observed in reality. Its movement 
can obey airflows or electromagnetic forces — it is of no importance. What means the 
explosion of fireballs? Some leakage of energy does exist, especially, by contacts with 
charged bodies. Balance of energy which is necessary to retain the current is failed, and 
plasma is extinguished. More exactly, the charge of the kernel becomes unsufncient to 
maintain electrons, they begin to fall into the kernel and this falling decreases still further 
its charge. The process becames avalanche-type and the fireball is quenched instantly. 
The pressure has no time to level off. The difference in pressure stimulate a shock wave 
like that of the spark or of the usual lightning. This wave is perceived as an explosion. 
The known phenomenon — the breakdown of the superconductivity by sufficiently strong 
magnetic field — explains nonstandard interactions of fireballs with electrical equipment 
[23]. 

The fireball is something like the model of the atom: a kernel and electrons spreaded 
on orbits. Of course, we use the term "geodesic flow" as a figurative expression. In 
accordance with the quantum theory and with the Heisenberg's indeterminacy principle, 
electrons are defined not as moving points but as a probability distribution. Nevertheless 
the term "geodesic flow" reflects much more adequately the bound state of electrons under 
the influence of powerful electrostatic field of the positively charged kernel. 

At latest works on the creation of fireballs [11] it was generated plasmoids (that was 
regarded by the authors as fireballs). It has been found experimentally that: "It was 
formed the negatively charged layer around the polarized kernel" . 

It may appear that the geodesic flow must lead to perpetual collisions of electrons. 
But the volume of the cover is macroscopic and it is huge by comparison with the size 
of an electron. Though any two big circles of S 2 intersect but electrons moving along 
these orbits (no matter how many they will be) can avoid collisions. Moreover, in the 
approximation when electrons are considered as points, the probability of collisions equals 
to zero. But the main reason: each separate electron cannot scatter by himself being 
bounded coherently with other electrons generating joint "plasma-sounded" vibrations. 
Thus the approximation of the geodesic flow is reasonable. 

Let us remark that all attempts to confine plasma was based on ring-shaped regions. 
But it needs strong magnetic fields. May be it would be better to try electrostatic fields. 
If there exists superconductivity of plasma, then it must help. Attempts of confinment 
are very difficult to realise because one endeavour to feed plasma by energy, to clamp it 
down, to make a pinch. However we may give freedom to plasma, to set it in a situation 
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of "singing", when plasma will have its own nondamping "sound" vibrations, moreover, 
vibrations which resist attempts of attenuation. Then it would be a confinement and a 
source of energy. It is exactly the state of superconductivity of fireballs. 

The Kapitza's conjecture on the nature of fireball was not exactly correct, because his 
experiments with SHF-vibrations gave a ring-shaped plasma rather than a ball. This was 
quite natural because of the lack of the kernel. But the idea of feeding is very important. 
One needs only firstly to create "atom", and after that to feed it by energy. 

One may imagine a mechanism of the fireball birth in the following way: suppose, for 
instant, that the curl of positively charged particles falls on the course of usual lightning. 
Then the negatively charged plasma may surround this curl, and the fireball being stable 
begins its life. This idea leaves room for experiments. For instance, to blow a ball of 
ionized, positively charged gas with the electron beam, or to lay it into the negatively 
charged medium. If it has been possible to create fireballs artificially then an incorporation 
in it exact doses of "fuel" and an energy- rejection would give a safe atomic reactor. 

So far as in the picture of superconductivity of Bardeen-Cooper-Schrieffer (BCS) the 
key part play phonons (quasi-particles related to the fundamental frequences of the cristal 
lattice), our goal is to find eigenfunctions of operator that describes oscillations of spherical 
cover of plasma into the centrally symmetrical electrostatic field and to prove the existence 
of a gap in the energy spectrum which provides superconductivity. 

Thus we take as a base the geodesic flow on the sphere S 2 , that is the uniform 
distribution of electrons on the bundle T^(S 2 ) of tangent circles of radius 1 to the two- 
dimensional sphere of radius R . The geodesic flow of electrons must be stable relative 
to small disturbances, since, due to Coulomb interaction, electrons resist compressions. 
Deviations from this stable movement provide oscillations which play the role of phonons. 
Interaction of phonons with electrons is the source of superconductivity. One needs to 
interpret elementary vibrations around the ground-state as independent oscillators and to 
replace it by operators. As in an atom, the potential of the kernel exceeds the interaction of 
electrons and it may be thought that electrons move along great circles with the constant 
speed. That was the reason that we begin with the bundle T^(S 2 ) of tangent circles to 
sphere. 

The kernel of the fireball consists of positively charged ions. The external cover of the 
fireball is the negatively charged plasma in which there is a current — the geodesic flow of 
electrons. That means that the density of electrons is uniformly distributed on the sphere; 
their velocities are uniformly distributed on the tangent directions to the sphere and its 
absolute value equal to the speed of circular motion around the kernel in response to its 
Coulomb attraction. As a result, the potential of the external cover acting upon internal 
points of the kernel equals to zero. This provide the stability of the kernel. Indeed, if an 
ion penetrates into the cover from the kernel then the electrostatic potential of the cover 
draw it back, since the potential of the external (relative to this ion) part of the cover 
equals to zero, meanwhile the potential of the internal part of the cover draws it to the 
center. Electrons located outside of the cover can not penetrate to the kernel, since the 
cover obstructs this penetration. Namely, though the total charge of the fireball is neutral, 
but in the vicinity of it, the negative charge of the cover is prevailing that repuls closely- 
spaced negative charges preventing its penetration into the kernel. Thus the positive 
charge of the kernel remains practically unchanged. So, to explore the kernel one can 
use the theory of nuclear matter [5]. Electrostatic potential sustains the movement of 
electrons in the cover which vaguely resembles the skin-effect in superconductors. 
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The fireball is something like soap-bubble. Soap-bubbles usually slowly fall down in 
the air, because the pressure inside it (caused by surface tension), and hence the density 
inside it, is a little more than outside. So the weight of the bubble is a little greater than 
the weight of the air. Similarly behave the fireballs. The additional commplication is 
that the kernel tends to spred due to Coulomb forces, and this electrostatic pressure is 
added to the usual pressure in gas. The resulting pressure turns out to be a little greater 
than atmospheric pressure which bring down the density in the central part of the kernel. 
Meanwhile in the boundary part of the kernel (at the vicinity of the cover) as well as in 
the cover itself the density must be greater due to the pressure in the central part. As a 
result the mean density of the fireball is only slightly greater than the density of the air. 
Consequently the fireball can freely float in the air with the slight tendency toward falling. 
The difference in pressure is responsible for the effect of surface tension as well as for the 
explosion of fireballs and the arising of shock waves when fireballs are extinguished. 

The main question remains to be answered. The question about losses of the energy 
by scattering. Take as a pattern the BCS-theory of superconductivity. In plasma of the 
cover there exist oscillations like that of Bloch-waves in cristallic lattice. The difference 
is that Bloch-waves are linked with the periodic potential of the lattice, while collective 
oscillations of the cover are caused by the spherical form of configuration space. The 
attracting field of the kernel and the difference in the pressure inside and outside of the 
cover make the cover resemble a spherical elastic film in which arise spherical oscillations. 
Further on, the mechanism of plasma stability becomes analogous to that of usual super- 
conductivity: scattering of electrons excite oscillations of the cover and that in its turn 
stimulate electrons. This process can be regarded (in the quantum level) as an exchanging 
interaction of electrons by using the analog of phonons. Our aim is to find the energetic 
spectrum of the secondary quantized problem and to detect a gap in this spectrum. 

Let us recall that the notion of current in our case is different from the usual elec- 
trodynamics notion. The movement of electrons is uniformly distributed on all tangent 
directions. So in place of a current one has only "microcurrents" which superimpose and 
do not give electromagnetic fields. Consequently, the energy does not release: the fireball 
does not radiate electromagnetic waves and there are no losses of energy for that process. 
The fireball only fluoresce as usual plasma. 

1. Electronic cover 

Consider at first vibrations of the density of the geodesic flow on a sphere. Imagine 
the cover as a homogeneous flux of a compressible fluid on the bundle T}(S 2 ) of tangent 
circles at the two-dimensional sphere of the radius R . On the sphere arise circular waves 
that can be parametrized by its centers — points of emittence of waves. At an initial 
moment, take the disturbance of the density being the S -function at a point O with 
velocities uniformly distributed on a tangent circle. As far as the velocity of superfluous 
density is always directed outside, the wave bend the sphere passing the equator, focusses 
to the diametrically opposite point (—0) turning again into the S -function. We will call 
such waves elementary. If a wave is not elementary then in place of focal points appear 
another caustics. Each elementary wave on the bundle T}(S 2 ) of tangent circles to the 
sphere is the torus Tq . Indeed in the initial moment one has the circle S* 1 at the fiber 
over the point O . Further, it is projected on a circle of the sphere S 2 with the unically 
defined velocities which is orthogonal to the circle. After two focusing, the initial circle is 
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identically maps to itself. Any two such torus intersect by a pair of anti-directed geodesies, 
that connect centers of the corresponded waves. 

Besides these waves of condensation there exist similar waves of rarefaction. To quan- 
tize means to correlate an operator with each elementary wave. Let us find the Hamilto- 
nian for corresponding periodic oscillations. 

Note that the potential energy of the flow (let us denote it by U) does not change 
during these oscillations. Indeed, from the point of view of geodesic flow, the disturbances 
at each moment are condenced at points of a circle which moves without deformation with 
the constant speed along the torus T 2 . The degenerations of the projection to a point in 
poles does not connect with the geodesic flow as itself, but with its projection on the base. 
More exactly, T^(S 2 ) — the bundle of unit tangent vectors to S 2 — is diffeomorphic to 
the real three-dimensional projective space 1RP 3 . Tangent vectors to the geodesic flow are, 
say, left-invariant vector field on MP 3 or on its two-fold covering S* 3 , that may be realized 
as the group of unit in absolute value quaternions. Under the stereographic projection 
from the center of the sphere S 3 on M 3 , trajectories of this field transfer into a family of 
Clifford parallel lines in Ell{3) which is the space M 3 equiped by elliptic metric induced 
from S* 3 by this projection. Shifts along the trajectories of this family are isometries. To 
any fiber S* 1 G T^(S 2 ) , i.e. to the unit circle in the tangent plane at a point I of the 
base S 2 , correspond quaternions that commute with I : the circle in S 3 passing through 
I and the neutral quaternion. It is a greate circle and so it projects onto a line in Ell(3) 
passing through the image of the neutral element. Under the Clifford shift, these lines 
run along a toric surface. Hence, the set of trajectories of points of a fiber over some 
point x G S 2 is two-dimensional torus. As was shown by Clifford, the metric of this 
torus induced by the ambient elliptic metric is Euclidian. Thus the potential energy of 
an elementary wave remains constant. 

Meanwhile, in contrast with the movement in a flat space, the momentum of the wave 
changes. This is caused by the fact that the movement is not free: it takes place not 
in the whole space but along the surface of a sphere. Take the spherical coordinate on 
the sphere with the pole O . Fronts of wave centered in O are circles ip = Const. As 
an momentum of this wave it should be taken the operator of differentiation by ip . But 
since this operator must be defined at all points (in particular at poles), it is necessary 
to multiply it by a function which is equal to zero at poles. The simplest such function is 
cosip. Thus as an momentum we consider the operator 

.h , d 

We will use the rational system of units in which the light velocity c and the Plank's 
constant H equal to 1. So, the Hamiltonian is defined by the operator 

Hf = -cosip^(cosip-^)f + Uf, (1) 

where U up to multiplicative constant corresponds to potential energy. We are interesting 
in solutions depending from ip only, so / = f(ip) ■ The equation for oscillations is 

d 2 f df 
— cos 2 ip-r-r: + cos ip sin ip— + Uf(ip) = 0, 
dip 2 dip 
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or 



L = 



d 2 f 



— tan^— — 



df U 



/ = 0. 



(2) 



# 2 



cfa/> COS 2 ^ 



We seek for the spectrum of the operator L defined by the equation (2). Eigenvalues 
of L will be denoted by a . We obtain the equation 



since the function / has to be independent on at poles. Let us note that the boundary 
value problem (3)-(4) defines a selfajoint operator, because the equation (3) can be written 
in the form 



Eigenfunctions of the selfajoint operator L\ with the boundary conditions (4) relative 
to the positively defined selfajoint operator of multiplication by the function (cos^) are 
mutually orthogonal on the interval [— n/2,ir/2] with the weight (cos^>). 

We found oscillations of density or the "longitudinal" oscillations on the sphere. Let 
us find its "cross" oscillations, that is the oscillations of the sphere itself from purely 
phenomenological point of view. With this in mind, let us solve the following problem 
that has an independent interest. It may be interpreted as a problem of oscillation of the 
balloon cover. 

Consider an elastic spherical film surrounding a volume with the internal pressure 
greater than that of the outer air. We seek for its eigen oscillations due tu surface tension. 

On any element of the surface act forces of surface tension which are applied to the 
boundary points of the element. The resultant of two forces which act on the opposite 
points of a given direction at a point O is normal to the surface and equals to the normal 
curvature of this direction. The resultant of all forces is the normal to the surface which 
equals to the integral of normal curvatures, hence, equals to the mean curvature at a 
point O . For the stationary state, these forces are balanced out by the difference in 
pressure. The difference between rersultant and the force caused by the difference in 
pressure (acting in the same direction) give the force acting on the cover. As far as we 
know, the problem in this statement was not considered. This model is different from 
usual statements of problems in elastic theory since usually one considers forces to be 
proportional to displacements. 

Let us formalize the problem obtained. Introduce spherical coordinates on the sphere 
of radius R . 



(cos 2 ip)f" — (sin ip cos tp) f' + (a cos 2 tp — U)f — 0. 
Boundary conditions for the equation (3) are 



(3) 




(4) 




{R cos ip cos 0, R cos ip sin 0, R sin tp}. 
Denote the corresponding point of the unit sphere by 



r = {cos -0 cos 0, cos ip sin 0, sin-0}. 
Denote a displacement of the point R,ip,<f> by 
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RU = {^,<p),R V ^,<p),RC(n}- 

The coordinates of the displaced point are X = Rr + RU . Let us introduce the 
moving frame on the unit sphere 

r = {cos -0 cos 0, cos -0 sin 0, sin if)}, 

A; = {-sin0,cos0,O} = ^, 

I = {— sin -0 cos 0, — sin -0 sin 0, cos-0} = r^. 

One has 

r<p<p = {cosip)^, k<f> = (-cosip)r + (sin0>)/, k^ = 0, ^ = -r; 
1$ = (-simp)k, r H = 1$= (sin^)/c, = = -r. 

Its vector products are 

[r x k\ — I, [k x I] = r, [I x r] = k. 

Calculate the mean curvature of the displaced surface, neglecting members of second 
order relative to II and its derivatives 

X$ = (Rcosip)k + RU^, X^ = Rl + RU^. 
The coefficients of the first quadratic form ds 2 equal 

(X^X^) = R 2 (cos 2 i) + 2 cos 7p(k, II )), 

(X^X^) = R 2 (cosip(k, n^) + (/,n^)), (5) 
(X^X^)=R 2 (l + 2(l,U ps , i )). 

[Xj x X$\ = R 2 {(cos ip)r + (cos^)[/c x U^} -[!x^]}. 
\[X^ x X^\ \ = R 2 ^cos 2 tj} + 2 cos 2 4>) + 2 cosijj(k, n ). 

1 - 1 {1 _ (1 , n ,)-*M}. 



| [X^ x X^,] I i? 2 COS "0 COS "0 

Normal to the surface is defined by the expression 

n = r + [jfc x n„] - [ 1±M - (/,n^)r - (6) 

COS-0 COS-0 

We have 

X H = R{(-cos 2 ip)r + (cos sin ^)/ + IL^,}, 
X H = R{(-smip)k + n 0V ,}, 
X^p = R{—r + IL^,}. 

The coefficients of the second quadratic form dN 2 equal 

{X H ,n) = R{-cos 2 ip- cos 0> sin 0>(r, IL^)) + (r,II^)}, 

(X H ,n) = i?{tan^(r, II ) + (r, IL^,)}, (7) 

(X^p, n) = R{-1 + (r, 11^)}. 
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Denote (2 x 2) -matrix of the quadratic form dN 2 — \ds 2 by A(X) . Using (5) and 
(7) find entries of this matrix 



A n = R{— cos 2 ip — cos ip sin ip(r, U^) + (r, 11^)} — Ai? 2 {cos 2 ip + 2 cosip(k, II^)}, 
Ai 2 = i?{tan ^(r, IT ) + (r, 11^)} - XR 2 {cosi;(k, 11^) + (/, n^)}, 
A 22 = R{-1 + (r, n^)} - Ai? 2 {l + 2(1, n^)}. 

Roots of the equation detA = define main curvatures. By dividing half of the 
coefficient of A in this equation by the first coefficient and neglecting members of the 
second order relative to II we find the mean curvature 

K = i (-i - ^<r, iw + £M + (liSal + (/ , n +) + *M . 

it [ 2 2 cos^ xp 2 cos ip J 

The force which act on the surface at a point (0, ■0) is directed along the normal 
n (see (6)). It is proportional to K with the coefficient of propotionality which we 
denote by (3 . The inner pressure gives the force which is proportional to j^n . To have 
an equilibrium for the case when the displacement of surface equls to zero (sphere), the 
coefficient of proportionality must be also equal to (5 . The summary force is equal 

^{-^( r ,n,) + ^4 + ^i + (i ,n,) + ^}. 

[ 2 2 cos 2 ip 2 cos-0 J 

One obtains the equation 



j£° = l3 r ( { -^V,n,) + p**L + M**) + (i ,n,) + . (8) 

ot [2 2cos 2 -0 2 cos-0 J 

Rewrite the equation (8) in the moving frame coordinates: u — (r, II), v — (k, II), w — 
('•")• 

1 d 2 « _ 1 f tain/) I m^ + (cos 2 VQm^ . I 1 

/3 <9t 2 R\ 2 2cos 2 V "T «A/> 1" cos y, J ■ 

19!" _ (J, &IL\ _ r, 

a t 2 - l^, 9t 2 ) — u, 
at 2 ~~ 9< 2 / ~~ u - 

If at the initial moment v — 0, v — and w = 0, u> = , then v = 0, w = . Thus, 
the linearization leads to purely radial vibrations. The final equation is 

1 d 2 u _ 1 [ tan?/; + (cos 2 jj)u^ \ 

(3dt 2 ~ R\ 2 2cos 2 ^ J " 1 j 

As expected, the right hand side of the equation (9) gives the Laplace-Beltrami oper- 
ator on the sphere. 

As soon as the equation (9) is autonomous relative to , we find the solution in the 
form of traveling wave relative to , that is u — u(4> — nt, ip) . We have 

k 2 if tan^ ua,6 , um\ 



(3 vv R{ 2 v ' 2cos 2 V 
Separation of variables u = <£>(0)\I/(-0) gives 
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= (- sin^ cos^)^ + (cos 2 )^ = 
$ 2Rk 2 cos 2 

The period of the function <£> must be equal to 2n , hence C = — n 2 , n e Z. The 
equation for \1/ has the form 

(cos 2 ^)^" - (sin ^ cos -0)*' + (2n 2 Rn 2 cos 2 ^ - n 2 /?)^ = 0. (10) 
The boundary conditions for the equation (10) are 

= *(-£) = <>, (ii) 



,2/ V 2, 

since the function \l> does not depend on at poles. The expression a = 2n 2 RK 2 can be 
considered as the spectral parameter for the problem (10)-(11). Zonal spherical harmonics 
are its solutions. 

The spectral problem (10)-(11) appears exactly the same as (3)-(4). The striking fact 
that entirely different approaches and different phenomena lead to the same equation 
(10) = (3) give confidence to the Hamiltonian (1). 

Eigenf unctions u m ^ of the Laplace-Beltrami operator on the two-dimensional sphere 
are enumerated by pairs of integers (m, I), \m\ < I and relate to the eigenvalues —1(1+1) . 

In accordance with the established linguistic tradition, we shall call particles which 
are obtained by quantization of the above-mensioned waves (both langitudinal and cross 
waves) by spherons. 

Theorem 1 The solution u = of the equation (9) is stable in Liapunov's sense. 

Proof. The equation (9) can be rewritten as the ordinary differential equation in the 
space / 2 

= < 12 > 

The spectrum of the operator A , which stand at the right hand side of the equation 
(12), equals 

A m ,z = -1(1 + 1), I > \m\. 
The equation (12) can be written in normal form as a system 

^ = v 

| = 1. < 13 > 



The spectrum of the matrix 



/ 

a o r 



standing at the right hand side of the system (13), is purely imaginary. It is equal 
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±iyjl(l + l). 

Indeed, if f k , k G N is a basis of eigenvectors of the operator A responding to 
eigenvalues — , then 

u = f k e ± ^\v = ±t\ k f k e ±lXkt 

is a fundamental system of solutions to the equation (13). 
The theorem is proved. 
□ 

2. Second quantization 

We will quantize the equation obtained in momentum space. The role of the momen- 
tum space will plays the space of spherical harmonics. Let us note that it has the discrete 
basis that may be defined by eigenfunctions of the Laplace-Beltrami operator. The thech- 
nique of quasi-discrete representation, i.e. the preliminary consideration of the problem in 
finite parallelepiped with periodic boundary conditions and after that the passage to the 
limit from discrete case to the continuous momentum representation [6] (one of the feeble 
item in the ground of quantum statistics), is found to be unnecessary due to compactness 
of the phase state. 

Denote the spherical harmonics by v mj i(s), s G S 2 . Expand a solution to the equation 

where A is the Laplace-Beltrami operator on the sphere S 2 , in the series in terms of 
spherical harmonics v m j(s) . Coefficients of this series are defined by the expressions 

u m ,i(t) = J u(s, t)v mi i(s) ds. (15) 

Apply the operator J^- to the (15) taking into account that the Laplace-Beltrami oper- 
ator is selfajoint and that the functions v mi i(s) are ^ s eigenfunctions with the eigenvalues 
-1(1 + 1) 



d 2 u m ,i(t) 



/ Au(s, t) v m>l (s)ds = / u(s, t) Av m< i(s)ds = -1(1 + l)u m j(t). 



dt 2 

One obtains a decomposed infinite system of ordinary differential equations 

^^ = -l(l + l)u m M (16) 
Each equation of the system (16) is the equation of harmonic oscillator with the unit 



mass and the frequence yl(l + 1) . We shall call it by the oscillator (m,l) . Apply to 
it the standard proceedure of quantization. That leads to second quantization of the 
oscillation of the geodesic flow of electrons. 
The Hamiltonian of the oscillator (m, I) is 
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Hl^\ip 2 + l{l + l)q 2 ). 

We enumerate eigenf unctions of the operator H^ t by indices n and denote it by 
£ n (m, I) (the number of an eigenfunction is denoted by the lower index). It is well known 

that the corresponding eigenvalues are E n = (n + \)\Jl{l + 1) . Introduce the operators 
[7] acting on eigenfunction as follows 



a + {m, l)£ n (m, I) = \Jn + 1 £ n+1 (m, I), a (m,l)£ n (m,l) = y/n£ n -i(m,l). (17) 

The ground-state of the system corresponds to n — , since in this case the next 
conditions are fullfiled a~(m,l)^o(m,l) = 0, and 

Urn,l)= {a+{r pl )r Um,l). (18) 
vn! 



The formulas (17) show that states £ n (m, I) are eigenfunctions of the operators a~ (m, l)a + {m, I) 
and a + {m,l)a~ (m,l) with the eigenvalues n+1 and n correspondingly. Consequently, 
it is natural to call the operator a + (m, l)a~(m, I) by the operator of number of spherons 
of type (m,l) . Denote this operator by h s {m,l). The commutation relations take the 
standard form 

[a~{m, I), a + (m, I)] = 1, [a~(m, I), a~(m, I)] = [a + (m, /), a + (m, I)] = 0. 

Meanwhile operators that relate to noncoinsiding values of indices commute. The 
Hamiltonian for an oscillator is H°(m, I) = \Jl{l + l)(a + (m, l)a~(m, I) + \\ s {m, I)) , or 

H°(m, I) = ^l(l + l)(n s (m, I) + h s (m, I)), 

where P(m, I) is the identity operator on the space of states of the oscillator (m, I) . The 
full Hamiltonian of the field of the cross waves is represented by the sum 

H°= £ H°(m,l)= J2 y/l(l + l)(n'(m,l) + ^l a (rn,l)). (19) 

l,\m\<l l,\m\<l 

One can express the Hamiltonian in the representation of occupation numbers. For 
brevity sake, let us relabel Fourier coefficients (m, I) of spherical harmonics by unique 
index h and denote the corresponding operators by a + (h) and a~(h) . The state of the 
system is characterized by a set of occupation numbers hi, hi, ...hN , i-e. the state 

l<k<N V^fc- 

means that there are hi particles of the sort . 

The mechanism of the superconductivity of plasma is somewhat different from that 
of electron-phonon interaction in BCS-theory. The later is connected with interactions 
of electrons and vibrations of ion lattice, that is with something external relative to the 
electron flow itself. Whereas electron-spheron interaction is the interaction of electrons 
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with vibrations of the flow of electrons. In other words, this is something like the selfcon- 
sistency of Hartree [26] . The only difference is that in Hartree theory the movement of an 
electron is consistent with its own field, wheareas in plasma, the movement of an electron 
which affect the creation of the field of the flow must be consistent with the field of the 
flow as a whole. In this situation the field of the flow is not the electromagnetic field, but 
the field of spheron's oscillations. Roughly speeking, fluctuations in state of an electron 
as if " damp" the electron flow giving to it quants of energy. In respons, the coherent flow 
slightly "pushes" this electron returning it back into "formation". 

As the components of the field of the flow it may be taken longitudinal and cross 
waves. The equation (3) describes longitudinal oscillations in the absence of cross ones; 
the equation (10) describes cross oscillations in the absence of longitudinal ones. But 
in the general case longitudinal and cross oscillations are mutually connected. Local 
deformations of the sphere lead to the changes of density; local decreasings of density 
decrease local forces of surface tension that leads to deformation of the surface. Note that 
the increasing of the function / decrease the function \P and vice versa. In other words 
oscillations of the longitudinal and cross waves are in opposite phase. Taking into account 
that these oscillations satisfy to the same equation, it will be supposed that \I/ = — / . 
So spherons will be described by a scalar field with the Hamiltonian (19). Note further, 
that an elementary wave of density arrive at the diametrically opposite point being in 
the same phase as at the initial point, that is the excess of density at the initial point 
leads to the excess of density after focussing. Thus we will be interested only in even 
functions relative to the antipodal mapping. Hence, the spectral parameter / must be an 
even number. The corresponding system of eigenfunctions are complete in the projective 
space MP 2 . The Hamiltonian of the field of spherons takes the form 

H°= \/2/(2Z + l)(n s (m, 21) + -I s (m, 2/)). (20) 

l,\m\<2l 

Our aim is to write the Hamiltonian which consists of 

1) . the free Hamiltonian of the spherons field (20), 

2) . the free Hamiltonian of the field of electrons, and 

3) . the interaction Hamiltonian. 

3. The gap 

We wish to reduce the situation to the case of the Bose statistics as in usual theory of 
superconductivity. In that theory it was used quasi-particles — Cooper's pairs. There is 
some mystery in the fact that a pair consists of electrons with the opposite momentum and 
spins. It is rather strange that particles flying in opposite directions are linked. Namely 
Cooper that proposes this idea calls this connection "deeply mysterious". Though, there 
was some statistical arguments in favour of it (see, for instant, [29]). The main reason is 
that the combined spin becomes zero which leads to the Bose statistics. So in the state 
with minimum energy there are plenty of particles, so called coherent condensate, having 
high-capacity inertia which obstacle small changes. Together with the energetic gap this 
gives the intuitive explanation of the superconductivity. 

In our case the combination in Cooper's pairs is much more natural than in BCS- 
theory. We again combine electrons with the opposite spin and momentum. But now 
electrons of the pair move along the same geodesic in the same direction. It is convenient 
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to think on it as situated at the diametrically opposite points. More exactly, the support 
of the function that describes the joint state of the pair belongs to the antipodal diagonal 
of the direct product S 2 x S 2 , that is it belongs to the manifold {x, —x} G S 2 x S 2 . In 
particular this means that if one of the electrons was detected at the point x then its 
pair will be at — x . Thus if one of the electrons of the pair jumps from the state with 
the momentum p into that with the momentum k then the wave of density having the 
concentration at diametrically opposite points transfers its pair from the state with the 
momentum — p into that with the momentum —k . 

Let us show that combining electrons in Cooper's pair provides the gap in energy 
spectrum. Denote the width of the electronic cover by r®. The phase space for electronic 
cover related to spheres with radii from R to (R + tq) is the direct product of the man- 
ifold of geodesies on a sphere (that is MP 2 ) by the interval [R, R + ro] . An momentum 
k and a spin, taking values ±1/2, parametry eigenf unctions. The set of eigenfunctions 
plays the part of a crystalline lattice being something like a framework maintaining super- 
conductivity. Accidental disturbance in density of electrons spreads along the sphere as a 
spheron and compensates the related disturbance. It is possible since the phase velocity 
of waves is far much than that of electrons. The kinetic energy of an electron is defined 
by the radius of the sphere on which it moves. Aside from interchanging processes due 
to spherons, the potential energy is defined by interelectron repulsion which we will often 
neglect in what follows. Interchanging processes give the effective attraction, and this at- 
traction is far-acting since oscillations of electronic cover is a global process acting on the 
whole sphere. So the interchanging processes give a negative contribution into potential 
energy of the system. That is the physical reasson for arising of the energetic gap in the 
spectrum. Let us note that the role of Debye wavelength which defines the characteristic 
linear length of the system, plays the width of the electronic cover r . 

Denote by Q the charge of the kernel of a fireball, and by R its radius. Let p be the 
average volume density of ions, and e its charge (in the case of single ionization). One 
has Q = ^nR 3 pe . Let p be the average volume density of electrons in the cover. The 
volume between spheres of the radii R and R + r equals 

4 4 
-ir(R + r) 3 - -irR 3 . 
3 v 1 3 

The attractive force to the kernel acting on an electron that lie at the distance R + r 
from the center of the kernel equals r^~p — t(R+r) 2 • ^ ut this f° rce i s shielded by the 
negative charge of the layer between radii R and R + r . Denote the shielding coefficient 
(depended on r only) by «(r) . Then the shielded attractive force of the kernel equals 
K ^ t{R+r) 2 ' Besides this force, on an electron acts the force of repulsion from electrons of 
the layer (R, R + r) only, since the resultant force acting from the external layer equals 
to zero. The resultant force acting from the internal layer equals to that of repulsion from 
the center induced by the overall charge of the layer (R, R + r) . Hence, the full force is 
equal to 

4irpR 3 e? _ Uit(R + rf AnR 3 \ 2 

K[j 'h(R + r y ~ \ 3{R + r y ~ 3{R + r y ) ^ e ■ 

This force equals to zero at r = r , consequently, the equation defining r is 

K(r )pR 3 = ((R + r ) 3 -R 3 )p. (21) 
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ro 
R 



The equation (21) is in a sense the condition of neutrality of the fireball. If one denote 
= 5 , then the equation for 5 is 



5 3 + 35 * + 35= <?)P. (22) 
H 

Find the velocity v of electron moving at the orbit of radius R + r . The equation of 
movement at the plane (x, y) are 

v v 

x — (R + r) cos — 1; y — (R + r) sin — 1. 

R+r R+r 

By differentiation one obtains the centrifugal acceleration 

I? „,2 



a = — 



m R + r' 
where m is the mass of the electron. Hence, 



v = 



'F(R + r) 



m 



The kinetic energy of an electron is 



e(R + r) = = ^0-^[ K (r)p& - »({R + r) 3 - R% 



It follows that 



eat) = 

Besides, in view of (22), e(R + r ) = 0. 

As in [28], imagine all admissible states in which electrons with momentum and spin 
((p, 1/2), (— p, — 1/2)) are combined in pairs. The scattering of a paire ((p, 1/2), (— p, — 1/2)) 
on a pair ((/c, 1/2), (—k, —1/2)) means the transition from the state ip in which the 
cell ((k, 1/2), (—k, —1/2)) is free and the cell ((p, 1/2), (— p, — 1/2)) is occupied, into 
the state x being in opposite condition. Denote by x 1 the probability that a cell 
((p, 1/2), (— p, — 1/2)) is occupied. Then the probability that a cell ((/c, 1/2), (—A;, —1/2)) 
is free equals y\ = 1— x\ . The probability of the state ip in which the cell ((k, 1/2), (—A;, —1/2)) 
is free and the cell ((p, 1/2), (— p, — 1/2)) is occupied is the product of probabilities 
P$ = Xp(l — x\) . The amplitude of the state ip is = x p yu ■ Analoguously the 
amplitude of the state \ m which the cell ((p, 1/2), (— p, — 1/2)) is free and the cell 
((k, 1/2), (—A;, —1/2)) is occupied equals P x = x^y v . Define the energy of an electron 
with the momentum p by e(p) , and denote the matrix element of the potential energy 
of interaction for the process in question by W P:k . The full energy of electronic cover 
(discarding the Coulomb forces of electronic repulsion) is 



,k%p%kypyk 

£2 £ (p)^ _ W(J2x p y p )(Y,x k y k ). (23) 

P p,k V P k 

Here W P:k is changed by — W for momentum related to tangent bundle to spheres of 
radii in the interval (R, R + r ) , and W p>k = for other momentum. The factor 2 of 
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the first sum arised from the fact that electrons in the states both (p, 1/2) and (p, —1/2) 
have the energy e(p) . 

The ground-state gives the minimum to the energy E . Let us find this minimum. 



Theorem 2 The following relation is fulfilled at the point of minimum of the energy E 

2W 



i = E 



^A 2 + £ 2 (/)' 



where 



Proof. 

Differentiate the expression for E relative to xf 



The notation used in the statement of the theorem: 



dE 



2e(l)-W(Y,*m) 



0. 



gives the possibility to rewrite the preceding equation in the form 



2e(0V^-^ = A(l-2^). 
By squaring we obtain the quadratic equation relative to xf 

A 2 



X, 



4(A 2 + e 2 (/)) 
1 ± 



0. 



£ 2 (0 



\ A 2 +5 2 (/) 



(24) 



(25) 



The sign minus should be taken in this expression, since otherwise the right hand side 
of the equation (25) would be negative. So 



1 + 



e\l) 



\ A 2 + e 2 (/) 



We shall show below that, as in the classical theory of superconductivity, the value A 
corresponds to the gap in the energy spectrum of the system. Substituting this expression 
into the formula (24), we obtain the analog of the famous gap equation: 



e 2 {l) 



1/2 



or 



E 



A 2 + e 2 (/) / 
2W 



r v /a 2 +5 2 (o' 



(26) 
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The theorem is proved. 
□ 



Let us introduce v(e) — the density of states of electrons with the energy e . Using 
this notion rewrite the gap equation in an integral form 



W 



I 



*W v{e)de 

(R+r ) 2 V / A 2 +£ 2 ' 



The oscillation of the function v{e) on the interval (R, R + r ) is insignificant. We 
shall denote v{e) simply by v and it can be carried out of the sign of the integral. We 
obtain 

de 



rs(H Hp 

/ J . (27) 

Je(R+r ) a/A 2 + S 2 V ' 



Wv Je(R+r ) VA 2 + E 2 ' 

The integral at the right hand side of (27) monotonically decreases from +00 to as 
A runs from to +00 . Consequently, a solution A to the equation (27) exists and is 
unique. Let us find the explicit expression for A by calculating the integral (27). 

2 _ Arsh £(g)_ Arsh £^±l2) =Arsh ^ 2 



Wv A A 3A 

Hence, 



, 2 2ne 2 pR 2 . , 

sh — - = -f . 28 
Wv 3A y ' 



It follows 



Let us show that A is actually the gap in the energy spectrum. Having this in mind, 
we evaluate the contribution £(/) from the adding of the combined pair of electrons with 
momentum ±Z into the energy of the ground-state (when all electrons are combined in 
pairs). From (23) we have 

= 2e{l)x 2 -2Wx iyi J2x k y k . 

k 

Indeed, the first summond in this formula is the kinetic energy of the pair (/, — /) , and 
the second is the contribution into the negative part of the energy of the ground-state 
arising from the possibility of interaction processes of the pair in question with all other 
pairs (k, —k) . The coefficient 2 reflects the fact that the pair (/, —I) in the sum (23) occurs 
twice: in summing by k , and in summing by p . Itroduce the notation E\ = ^JA 2 + e 2 {l). 
Taking into account values: xf = \{\ — ^A) , and denotation A = WY,k x kUk , we find 



= 25(/)|(l- f )-2A 



2,» \V2 



e 2 (l) 
E'f 



*(0 " * - Af = e(l) - E h 
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If to add to the ground-state merely one (unpaired) electron with the momentum I 
then the pairs (I, —I) cannot contribute into the energy of the ground-state, and hence 
the energy of the system after such adding become 



E - f (Z) + e(l) = E + Ei = E + ^A 2 + e\l). 

Consequently, the adding of an isolated electron raises the energy of the ground-state 
at least on A . This means that to change the ground-state, i.e. to disjoin a combined 
pair, it is required the energy non less than 2A . This is the same as saying that there is 
a gap in the energy spectrum. 

To estimate the shielding effect n(r) , we will use the methodology of Tomas-Fermi, 
developed for computations of heavy metals spectra. The Tomas-Fermi model for centrally 
symmetric fields is based on the supposition that there are greate many electrons in a 
domain where the potential changes slowly, and in this domain the Fermi statistics is 
applicable. This supposition is fulfilled in our model of fireballs. 

Firstly, let us remined the deduction of the Tomas-Fermi equation in suitable deno- 
tations. Denote by po(r) the maximal absolute value of momentum of electrons situated 
at the distance r from the center. Let us agree that is the excess of the energy of 

an electron over the level of the termal potential. We have 

V(r) + ^ = 0. (30) 

Let an electron moves freely in a domain Q with the electron density jx = |y . Then 
the number of quantum states with the absolute values of momentum from p to p + dp 
is 

where Kk = p , and the factor 2 takes into account two possible orientations of the spin. 
We integrate from to k and equate the result to the full number of electrons N : 

or 

kl = 3ttV 

Let now suppose that a subdomain O C Q is selected being sufficiently large for 
correctness of (30) and sufficiently small, such that the potential energy in it changes not 
too much. Then 

By substitution (31) into the Poisson equation for electrostatic potential — V/e with 
the density of charge — e/i : 

AV = -47re 2 /i, 
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one obtaines 

= -3^( 2m)3/2 (-"') 3/2 

Using dimensionless variables 

r 2 7/3 rV 

X = T = To \9rt r i / = 



6 (3tt) 2 / 3 ' J Qe 2 ' 
where Q is the total charge of the kernel, one obtains the Tomas-Fermi equation 

dx 2 -Jx 



(32) 



Let us apply the equation (32) to our theory of fireballs. While calculating heavy 
metal spectra, the following boundary values for the equation (32) is customarily used 
(cf. [28], [4]) 

/(0) = 1; Bm/(x)=0. (33) 

However, we need for our model another boundary conditions. It must be taken that 
at the distant R we have the Coulomb potential, and the cover ends at points where 
termal fluctuations of electrons majorize the attractive forces. At the boundary points, 
the potential will be equal to the fixed value — ( , where ( corresponds to the termal 
energy of electrons at the temperature T . The related boundary conditions are 

f(R) = l; f(R + r ) = ^-^(. (34) 

Let us explore the geometrical structure of solutions to the equation (32). 

Theorem 3 The equation (32) admits the group G of scale symmetries allowing to re- 
duce it to a two-dimensional system which does not depend on the argument x . 

Proof. 

Let us rearrange the equation (32) as the system 

tjflll (35) 

dx \fx 

The group G can be found in the form / i— > A/, g i— > fig, inra. These transfor- 
mations transfer the system (35) into the system 

= vg 

v dx 

nil = ftff 3 ' 2 

V dx V V yfx 

To hold the initial shape of the system (35) it is necessary that A = fiu, fi 2 = A'V 
or, equivalentely, A = v~ z , fi = u~ A . Hence, the system (35) turns into itself if / i— > 
V ~ Z J ■, 9 l— v ~^9i x l— vx . It allows to reduce the dimension of the system. To do 
this, one needs to take as new variables such combinations of old ones that are invariant 
under the action of G . The simplest such combinations are y = x 3 f, z = x A g. In 
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order that the system in the new coordinates becames autonomous, (i.e. does not content 
the independent variable) we put x = exp p . In this case shifts of p transfers into 
multiplications by a constants v . Then the system (35) takes the form 

T- = 3y + z 

The theorem is proved. 
□ 

Now we can explore the phase portrait of the system (36) on the two-dimensional 
phase plane (y, z) . Besides the obvious equilibrium at the origin we have the singular 
point A = (y = 144, z = 3 • 144). By returning to the initial variables, the singular 
point A gives the famous partial solution to the equation (32) guessed by Sommerfeld, 
/ = 144a; -3 . It is agreed (cf. [4]) that this solution gives adequate asymptotics of the 
exact solution for atoms as r — > +oo , but does not satisfy the boundary conditions at the 
origin. This fact is understandable if to use the system (36). It is easy to check that the 
origin is the unstable knot of the system (36), and the point A is the saddle. The exact 
solution for the atom that satisfies both boundary conditions is the separatrix passing 
from the origin and tending to A as p — > +oo . It is not hard to show the existence of 
the separatrix. And for it r — > as p — > — oo . 

However we are interesting in boundary conditions (34). In terms of the system (36) 
these boundary conditions transfer into the following ones (let us recall that y = x 3 f, z = 

x 4 g.) 

2 7 R 3 2\R + r )\ 



where 



, 2 7 / 3 R , 2 7 / 3 (i? + r ) 

po = In pi = hi 



( 37r )2/3> ^ ( 37r )2/3 

Finally we have 

7,nn 2VBR \ — 27rB 

My 111 (3tt) 2 /3 ) - (3tt)2 ' 

My 111 (37r) 2 / 3 > — Qe 2 (37r) 2 S 



Thus we need a trajectory of the phase plane (y, z) that starts with the vertical line 
y = -j^fp- and leds to the other vertical line y = 2 q^?^°X C during the time p\ — p = 
lnr + | In . Coordinates of the final point (y(pi), z{pi)) of that trajectory give the 
derivative of the potential ^(R + r ) , and that in its turn is equal to the shielded force 
of attraction of the kernel. Because of z — gx 4 , g = jf, f = —jT2 , we obtain 



ftr' J Qe 2 

7 (n \ - (R+ r o) 4 9(rV) Or 
Z \I J D — b 4 9r 9x 



(37) 



Substitute into the equation (37) the value re(ro) from the equation (21). We obtain 
the relationship for finding the unknown parameter tq . Let us remark that it is natural 
to take for R values of fireballs that have been surveyed R ~ 10cm (cf. [23]). 
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The system (37), (21) can be solved by Newton's method. It is possible to take for ini- 
tial data nonshielded potential Vi(r) = —Qe 2 /(R + r); «i(r ) = 1 , and sufficiently small 
value r , say (r )i = -R/100. Find the corresponding value (po), (pi), 2/(Po)? 2/(Pi) • Us- 
ing (36) and reversing the time current, we obtain the trajectory yi(-), z\{-) . To calculate 
the derivative of the first approximation, we have to consider the variational equation 
of the system (36) for the trajectory yi{-),z\{-) . We obtain the derivative of K\ at the 
point /ti(r ) . Using Newton's algorithm, find k 2 and new (r ) 2 . After that the process 
is repeated. 

4. The field of electrons. 

Our aim is to quantize the field of electrons itself that moves at trajectories of the 
geodesic flow on the sphere. To simplify the construction we will take the tangent bundle 
of unit tangent vectors T l S 2 to the two-dimensional sphere as the phase space. This gives 
us the possibility disregard the electromagnetic field of the kernel, since its attraction was 
properly accounted. Objects of quantization are pairs of electrons, i.e. bosons having 
integral-valued spin. 

One of the way is to quantize the manifold T 1 ^ 2 , using the technique of geometri- 
cal quantization in terms of Kirilov-Kostant-Souriau [2], [14], [15], [16], [19]. [21], [22], 
[25], [27]. It would be very interesting to try this approach and to develop the tech- 
nique of geometrical quantization, considering bundles over the Kepler manifold and the 
projectivization of its sections as states. 

But we will use another approach. As far as electrons are moving at trajectories of the 
geodesic flow, we take as the main object of the theory (quasi-particles) geodesies itself. 
Consider the manifold of nonoriented geodesies on the sphere S 2 [3]. This manifold 
is MP 2 . Let us recall that the tangent bundle of unit tangent vectors T l S 2 to the 
two-dimensional sphere is RP 3 . The stereographic projection of MP 3 from the center of 
hemisphere induces on R 3 the metric of the elliptic space Ell(3) . The geodesic flow gives 
by this projection the set of parallel (in the sense of Clifford) lines. This set constitute 
the manifold RP 2 . The distant between Clifford-parallels remains constant in the metric 
of Ell (3) . It is natural to take it as a distant between geodesies on S* 2 . So, our quasi- 
particles are points of RP 2 . 

Let us find interrelation of quasi-particles. Denote the corresponding potential by 
V(p(x, y)) , where p(x,y) is the distant from x to y in RP 2 . To be uniquelly defined 
on RP 2 x RP 2 , the function V(p) must be zero at infinity, i.e. V(n/2) = , and besides, 
V(p) — > oo as p — > . This defines a dynamical system on a set of points of RP 2 . For 
finite sets of points it is the system of ordinary differential equations. Equidistant systems 
of points on RP 2 are stable equilibriums. These equilibriums are not isolated, since the 
set of points can move along RP 2 as a rigid body. If equidistant systems do not exist 
then we have no equilibriums. But, in view of condition V(p) — > oo as p — > , for any 
fixed value of energy, distants between points during the motion admits lower estimation 
by a positive constant. 

Thus, the phase space is the projective plane RP 2 (recall that spherons are parametrized 
by points of RP 2 also). To find the law of interaction of quasi-particles, let us calculate 
the interaction of two uniformly charged geodesies on S* 2 . 
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Theorem 4 The potential of interaction of two uniformly charged geodesies on S 2 equals 

1 



sin#' 

where 6 is the angle between these geodesies. 
Proof. 

First of all, find the interaction of a point A e S 2 with the uniformly charged geodesic 
<7i G S 2 . Introduce in S 2 the spherical coordinates x = cos cos ip,y = sin cos -0, -2 = 
simjj. Without loss of generality, take o\ as 

x = cost, y = sint, z — 0, (0 < £ < 27r). 

Take a point A with the coordinates (0, -0) . The vector I that joins A with the 
moving point t of a± , is 

/ = (cos cos -0 — cos t, sin cos -0 — sin t, sin -0) , 
|/| = ^2 — 2 cost cos cos -0 — 2 sint sin 0cos-0- 

The resultant force for all points of the o~\ is aligned with the meridian passing 
through the point A that is orthogonal to o~\ , i.e. along the tangent to the circle 
(x = cos cos s, y = sin cos s, z = sins at the point s = ip. This is the vector 
e= (— cos0sin-0, — sin0sin-0, cos-0). The angle between I and e will be denoted by 7. 
Then 

// \ cos t cos </> sin i/i+sin t sin </> sin ^> 

COS f — I L , c ) — > — 

\/2— 2 cos t cos flicos i/>— 2 sin t sin </> cos tp 

sinipcos(t—<p) 
^2—2 cos tp cos (t—<j>) 

The element of the force acting on A from the portion rft of the circle <j\ , equals 

, Idt Idt 
^' = W = (2-2cos^cos(t-0)) 3 / 2 ' 
The projection on the direction of the resultant is 

(p JF(tW = sin^cos(t-0) 

1 ' {)) (2-2cos^cos(t-0)) 2 " 

The force acting to the point A from all circle o~\ equals 

r-2-ir sin ijj cos(t— <j>)dt 
JO (2-2 cos ip COS (t— </>)) 2 

r7r sin i/> cos(-r)dT 
Z JO (2-2cosi/>cos(t)) 2 " 

Rearrange this integral by the change of variable tan \—u 



sin 



— "Co (1— u 2 )du 



2cos3 4 J , „ , , 

^ Sin ? j 1 1 \ _ 7T COS V> 

8cos 3| \ten^ tln~f J ~ 2sin 2 V' 



Find the interaction of the geodesic U\ with that a 2 obtained from it by rotation 
round the axis Ox by an angle 6 . 
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o" 2 = (cos r, sin r cos 9, sin r sin 9) . 
The spherical coordinates of a moving point r G cr 2 are 



tan0 T = cos6>tanr; cos0 T = —, = ; 

y 1+cos 2 Stan 2 t 

sin^v = sinr sin#; cos ip T = Vl — sin 2 r sin 2 9. 



sin0 T = 



cos 8 tan r 



cos 2 9 tan 2 r 



Let us remark that vl — sin 2 r sin 2 = V cos 2 r + sin 2 r cos 2 9 . The force F T with 
the absolute value 



it coster cos 2 r + sin 2 r cos 2 

r 2 sin 2 2 sin 2 r sin 2 

acts (along the meridian) on a point r We need its moment relative to the axis O x 
Hence, we project F T on the normal to the plane a 2 , i.e. on the vector 





f 1 





^ 




f ^ 


V = 





cos 6* 


— sin 9 j 


(!)■ 


— sin 9 i 




Vo 


sin 6* 


cos 9 J 




\ cos# / 



The unit vector £ tangent to the meridian passing through the point r , equals 

£ = (— cos0 T sin-?/v, — sin <p T sin^v, cost/v) = 

V cos 2 r + sin 2 r cos 2 9 ) . 



sm r sin f cos r 



sm t sin f cos ( 



\J cos 2 r+sin 2 r cos 2 9 ' ^/ cos 2 r+sin 2 r cos 2 ' 



Further on 



, sin rsin 0cosr + cos r^COS 2 T + SU1 2 T COS 2 # = 
y cos 2 r+sin 2 r cos 2 9 



\J cos 2 r+sin 2 r cos 2 9 

This means that we have the projection 



7TCOs9 



sin 2 r sin 2 ' 

The arm of this force is the distant from the point (cos r, sin r cos 9, sin r sin 9) to the 
axis O x . It equals sin 2 r . The moment acting on the portion dr of the circle 02 equals 

, , 7T COS 

M T dr = tt—ut 

2 sin 2 

The total moment acting on the 02 is 

r 27r , , , 7T 2 COS0 

M T rfr = — -2— 
sin 



r 

io 



Thus, points of the projective plane MP 2 are repelling (up to a multiplicative constant) 
with the force 

cos 9 
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where 6 is the angular distant between them. The potential of this force is 

1 



sin# 

The theorem is proved 
□ 

It is reasonable to believe that the kinetic energy of quasi-particles, as well as for 
spherons, is described by the Laplace-Beltrami operator. In the homogenous coordinats 
(x, y, z) on the projective space, this operator has the form 

D = (y 2 + z 2 )£ + (z 2 + x 2 )^ + (x 2 + y 2 )^~ 

2xy 1 g r - 2yz 7 g r - 2zx^- (38) 

a dxdy a oyoz dzdx v ' 

The distant between points of MP 2 is defined by the logarithm of the cross ratio. In 
view of the preceeding theorem, the potential of the pairwise interrelation of points with 
the coordinates {x^y^Zj) (xj,yj,Zj) is the sine in power —1 of the corresponding angle 



\j (a? + yf + z 2 )(x 2 + y] + z]) 



Vij = 9 I j 

sj {xiyj - x jyi ) 2 + (y iZj - y jZi ) 2 + (ziXj - z jXi ) 2 

where g is the interaction parameter. 

Denote by Di the Laplace-Beltrami operator in Ell(JBL 2 ) applying to the coordinates 
qi . The Hamiltonian of the system that consists of N interacting quasi-particles is 

N l 

^ 1 = -7E^A + E^- (39) 

i=l Z i<j 

As it was expected, the Hamiltonian is homogeneous. We consider it on the space 
MP 2 , i.e. on functions of zero degree homogeneity relative to Xi,yi, z^ for any i . Besides, 
it is invariant relative to the group of permutations of points. We obtain an analogue of 
the Calogero-Sutherland-Moser problem (CSM) [24]. But here the part of the circle plays 
the projective space MP 2 . It would be good to find the asymptotics of the spectrum as 
N -> oo . 

We will use the Euler formulas of the first and second orders for functions of zero 
degree homogeneity. 

x 9£ _i_ y d£ i z 9£ — n. 

X dx + y dy + Z dz U ' 



x dx? ^ y w w 1 ^ x y dxdy ^ y^dydz ^ zx dzdx u - 

By adding to and subtracting from D the value 



Q2 * 2 r,2 

2 u i 2 u I 2 

x + y Tr^ + z 



dx 2 dy 2 dz 2 ' 
and then by using the Euler formulas, we reduce (38) to the form 

q% \dx 2 dy 2 dz 2 ') 
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Hence, the operator H 1 takes the form 



-7 N 

^ = -^EI^I 2 ^ + E^'- ( 4 o) 



2^ 



Enumerate eigenvalues of the Hamiltonian (40) by natural numbers k and denote 
them by — a k . The corresponding complete orthonormal system of eigenf unctions is 
denoted by Y k (S), S G (S 2 ) N . Expand a solution of the equation 

w = H ^ u) (41 » 

in term of functions Y k (S) . Coefficients of the series are defined by expression 

£4(0= fu(S,t)Y k (S)dS. (42) 

J(S ) 

Apply the operator J^- to (42). Since the operator H 1 is selfajoint and Y k (S) are 
its eigenfunctions with the eigenvalues — a k , we obtain 



d 2 U k (t) 



[ H 1 U(S,t)Y k (S)dS = [ U(S,t)H 1 Y k (S)dS = -a k U k (t). 



dt 2 J(s 2 ) N ' ' ' J(s 2 y 
Hence, we get the disjoin system of ordinary differential equations 

dt2 = -VkUkit). (43) 

Each equation of the system (43) is the equation of the harmonic oscillator with the 
unit mass and the frequency x /a k (we shall call it by the oscillator k ) . The application 
of the standard quantization proceedure leads to the secondary quantization of the quasi- 
partical oscillations. The Hamiltonian of the oscillator k is 

We anumerate eigenfunctions of the operator Hi by indices n and denote by H n (fc) . 
It is known that the corresponding eigenvalues are 

E n (k) = (n + 7})V°k- 
We introduce operators as in (17). 



A + (k)E n (k) = V^+TZ n+1 (k), A-(k)E n (k) = y/nZn-^k). (44) 
The ground-state of the system corresponds to n = since 

A~(k)~ (k)=0, 

and 

Z n (k) = { -^^E (k). (45) 
vn! 
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The formulas (44) show that states S n (/c) are eigenfunctions of the operators A~(k)A + (k) 
and A + (k)A~(k) with the eigenvalues n + 1 and n correspondingly. Consequently, it 
is natural to call the operator A + {k)A~{k) by the operator of number of quasi-particals 
of type k . Denote this operator by n q (k). The commutation relations take the standard 
form 

[A-(k),A + (k)] = 1, [A-(k),A-(k)] = [A + (k),A + (k)]=0. 

Meanwhile operators that relate to noncoinciding values of indices commute. The 
Hamiltonian for an oscillator is 

H l (k) = ^S<Jk(A + (k)A~(k) + 

or 

H\k) = yfr(n' + ±V(k)), 

where l q (k) is the identity operator on the space of states of the oscillator k . The full 
Hamiltonian of the field of the quasi-papticles is represented by the sum 

H 1 = J2 H\k) = £ V^(n 9 + ^ q (k)). (46) 
k k z 



Let us find the operator H int of interaction of spherons with quasi-particles. The 
interaction annihilates a quasi-particle of an energy E 3 that corresponds to the oscillator 
k 3 , and this excites a wave of an energy E\ that corresponds to the spheron (mi, 2/i) . 
This wave create quasi-particle of an energy E 2 that corresponds to the oscillator fc 2 , 
and this annihilates a wave of an energy £4 that corresponds to the spheron (m.4, 2/4) . 
So, we can write the operator of interaction H mt . 



jjint 

-^E^u^Ja+K,^)^^^ (47) 

where are energies of the corresponding spherons and quasi-particles. The 5 -function 
at the right hand side of (47) provides energy conservation. 
The full Hamiltonian is 



H = H° + H 1 + H int = 
Ei,\ m \<2i pl{2l + l)(n s (m, 21) + ip(m, 2/)) + ^(n« + \l q {k))- (48) 
^E fe3 ,/ l ,| mi |< 2 ; 1 (a + (m 1 ,2/ 1 )A + (A; 2 )A-(A;3)a-(m 4 ,2/4))(5(E 1 + E 2 - E 3 - E 4 ). 



The problem is significantly simplified if we neglect terms of Coulomb interaction 
considering free quasi-particles as a first approximation. 

At first, take the special case when 7 = 1. In this case both Hamiltonians H° and 
H 1 act on isomorphic spaces that we denote by 21 s and 2l 9 correspondingly. Both these 
operators are diagonal in basis consisted of spherical harmonics. Its spectra coincide with 
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N-multiple spectrum of the harmonic oscillator. Note that operators and A ± com- 
mute, acting on different objects (different components of a state- vector) The Hamiltonian 
equals 



H = E*, H <2* yj2l(2l+l)(n'(m, 21) + §F(m, 2/)) + n«(m, 2/) + ±F(m, 2/))- (4g) 
W E^i^i^^a+K, 2Zi)A+(m2, 2/ 2 )A-(m 3 , 2l 3 )a-(m 4 , 2l 4 ))5{k + l 2 - l 3 - l 4 ). 

We shall call the interaction switch-back if mi = 777,3, ^1 = ^3 and m 2 = WI4, l 2 — l 4 . 
Find the spectrum of Hamiltonian under the natural conjecture of switch-back interaction. 
In this case the interaction behaves as if it were two different acts. The first act is the 
annihilation of a quasi-partical of type (mi, 2^) with the synchronous excitation of a 
spheron of the same type. The second one is the annihilation of a spheron of type (m 2 , 2Z 2 ) 
with the birth of a quasi-partical of the same type. Since the Hamiltonian contains the 
sum of all these processes, one can imagine that the state-vector of quasi-particles turns 
into the state-vector of spherons and conversely. The interaction Hamiltonian is 

H = Ei,|m|<2j ^21(21 + l)(n s (m, 21) + |P(m, 2/)) + 

E*, H < 2 * y/2l(2l +l)(n«(m, 21) + ±P(m, 2/))- (50) 
W E^K|< 2il (a + (mi, 2h)A + (m 2 , 2l 2 )A-(m u 2Z 1 ) a -(m 2 , 2l 2 )). 

Let us rewrite the operator H given by the formula (50), by using the polarization 
of the state-space 21 = 21 s © 21 9 . Take in both isomorphic subspaces 21 s and 21 9 the 
concordant bases. Then the operator H mt interchanges vectors of these subspaces and 
multiply that by — W , i.e. H mt : (x,y) 1— > (— Wy, — Wx) . Consequently, the operator 
H has the block form: 

( Z Z ) ■ («) 

where B u = B 22 = Diag(bi, b 2 , . . .), B 12 = B 21 = -WI . 

To find the spectrum of the operator H , we explore a little more general problem. 
Consider an operator A , acting on the direct sum of two isomorphic separable Hilbert 
spaces K © K , with the block structure 



A= Z\ , (52) 




where Bi are commuting selfajoint operators acting on the space K . In view of the von- 
Neumann theorem [18], there exists selfajoint operator C such that both operators Bi are 
functions of C . Take E\ — the spectral decomposition of the space K , corresponding 
to the operator C . Then 



£1 = J a(X)dE x , B 2 = J (3(X)dE, 



Note that the diagonal subspaces K\ = {(x,x)} and K 2 = {(x, — x)} are invariant 
relative to the operator A . By using the polarization K © K = K\ © K 2 , we reduce the 
operator A to the block-diagonal form 
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A = 



( 



(B 1 + B 2 ) 






(Si - £ 2 ) 



) 



(53) 



Denote by £\ © E\ the operators that correspond to the spectral family E\ into the 
subspace K 1 = {(x,x)} , and by E\ Q E\ the operators that correspond to the spectral 
family E\ into the subspace K 2 = {(x, —x)} . 

We proved the following theorem 

Theorem 5 The spectral representation of the operator A on the space K 1 is 



Corollary 1 The spectrum of the operator H can be obtained from the spectrum of the 
operator Bn = B 22 by the shifts (as a whole) both to the right and to the left for a 
distance of W . 

Now consider the case when 7 = n/m any rational number. As before we neglect the 
interactions defined by terms Vij . Let us unify n exciting states of spherons in a cluster 
that we call pseudo-spheron. In the same way, m exciting states of quasi-particles will be 
a cluster called pseudo-particle. For sufficiently many quasi-particles and quasi-spherons, 
the Hamiltonians H° and H 1 will asymptotically coincide. Irrational values of 7 can 
be approximated by rational ones preserving the asymptotics of the spectrum. So, the 
problem can be reduced to the special case when 7 = 1. 

Since the spectrum of the operator Bn = B 22 is known, then at the approximation 
in question we obtain the spectrum of the Hamiltonian H . So, we obtain a possibility 
(in principle) to find mean values of important operators. Besides, our model Hamil- 
tonian provides the basis for application of the perturbation theory for more realistic 
Hamiltonians. 

The question: 

Does there exist a coherence of electrons under a line (usual) lightning? To resolute 
the question on the quantum character of the usual lightning (or any spark discharge) 
one needs to explore oscillations of cylindrical flow of electrons. Propose some conjectures 
concerning this question. 

5. Conjectures 

The inner canal of a lightning is the ionized, positively charged plasma. Electrons 
are winding round this canal into a cylindrical cover. Stepped leader is, apparently, a 




The spectral representation of the operator A on the space K 2 is 




27 



process of transformation into a superconductive state. And the superconductive state 
itself is connected with macroscopic separation of positive and negative plasma. It seems 
that Abrikosov's threads in the theory of usual (low-temperature) superconductivity have 
the same structure: canal from ions surrounded by a cover of winding electrons (not for 
nothing these threads are called vortical). It is precisely this separation of plasma which 
gives rise to superconductivity. The ionized air in the course of thunderstorms behaves 
like a superconductor of the second kind. Stepped leader is the process of generation of 
analogs of Abrikosov's threads — media of plasma superconductivity. These threads can 
curve, branch, and terminate due to rushes of the wind. The lightning flying through these 
canals repeates its form. The stepped leader constructs the canal of lightning by steps, as 
if it makes links of a chain. So, one of the possible ways of forming fireballs is separation 
of one of the chain from the canal of lightning. Nice recent theory of lightning (Runaway 
Breakdown) [9] takes into account the leading role of high-energy cosmic rays. It seems 
that Runaway Breakdown relates rather to the stepped leader then to the lightning itself. 
This is in agreement with the recent observations [10], where it was shown that steps 
of the stepped leader were well correlated with the high-energetic microsecond burst of 
gamma emission. 

The question arises: why the lightning extincts so fast if there exists the canal of 
superconductivity? Perhaps it happens due to the magnetic field. Too strong magnetic 
field, that is excited by the powerfull current in the superconductor, destruct its creator 
- the state of superconductivity that induces this current. Electrons begin to fall from 
the cover into the canal, superconductivity disappears, and the lightning extincts. The 
difference in pressure generates the shock wave — the thunder. 
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